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We investigate the prospects of creating broad rotational wave packets by means of molecular 
interaction with long sequences of intense femtosecond pulses. Using state-resolved rotational Raman 
spectroscopy of oxygen, subject to a sequence of more than 20 laser pulses with peak intensities 
exceeding 10^® W/cm^ per pulse, we show that the centrifugal distortion is the main obstacle on 
the way to reaching high rotational states. We demonstrate that the timing of the pulses can be 
optimized to partially mitigate the centrifugal limit. The cumulative effect of a long pulse sequence 
results in high degree of rotational coherence, which is shown to cause an efficient spectral broadening 
of probe light via cascaded Raman transitions. 

PACS numbers: 33.80.-b, 05.45.Mt, 42.50.Hz 


I. INTRODUCTION 

Control of molecular rotation with ultrashort laser 
pulses is an active research area of experimental and the¬ 
oretical molecular science. New regimes of molecular dy¬ 
namics and molecular interactions can be accessed and 
studied by means of controlled rotational excitation, e.g. 
by varying the frequency of molecular rotation and con¬ 
trolling its directionality. The ability to create and shape 
rotational wave packets is instrumental in achieving high 
degrees of molecular alignment and, hence, in exploring 
the electronic structure of molecules. For recent reviews 
of this broad topic, see m- 

Significant progress has been made recently in extend¬ 
ing the reach of rotational excitation. An “optical cen¬ 
trifuge” has been used to spin molecules up to the 

extreme frequencies on the order of 10 THz [5] and to 
create coherent superpositions of more than 50 rotational 
quantum states [3- Owing to the adiabatic mechanism 
of the centrifuge spinning, controlling the phases of the 
individual states in such ultra-broad wave packets is not 
a simple task. An alternative route to generating ro¬ 
tational wave packets with well defined relative phases 
employs an impulsive excitation by an ultrashort laser 
pulse. 

The width of the wave packet created by a femtosecond 
laser pulse, acting as an instantaneous rotational “kick”, 
can be characterized by a dimensionless kick strength 
P = Aa/{4:h) J £'^{t)dt, where Aa is the polarizability 
anisotropy of the molecule, £ is the temporal envelope 
of the pulse and h is the reduced Planck’s constant [8]. 
Being proportional to the intensity of the laser field, P is 
limited by the ionization threshold of a given molecule. 
In a typical case (e.g. N 2 or O 2 ) this means that a wave 
packet of only a few rotational states can be produced by 
a single kick. 

The ionization limit can be avoided by breaking a sin¬ 
gle laser pulse into a series of N pulses (or a “pulse 
train”). If the pulses in the train are separated by the 
so-called quantum revival time Tj-ev, the cumulative kick 


strength P/v of the whole train is the sum of the kick 
strengths of its individual pulses [H]- Hence, broadening 
the wave packet can be achieved by increasing the num¬ 
ber of pulses, while keeping their intensity below the ion¬ 
ization threshold. Multiple schemes of using long pulse 
trains (i.e. A ^ 1) for the controlled rotational excita¬ 
tion of molecules have been proposed theoretically [SHI3], 
and implemented experimentally with sequences of up to 
eight laser pulses [HHH]. 

The revival time is inversely proportional to the second 
derivative of the rotational energy E{J) with respect to 
the angular momentum J [isiiini- In the case of a rigid 
rotor, Trev reduces to l/2cB (where B is the molecular 
rotational constant and c is the speed of light in vacuum) 
and does not depend on J. However, fast molecular rota¬ 
tion leads to a non-negligible centrifugal distortion and a 
J-dependent revival time. As a result, the accumulation 
of the total kick strength from pulse to pulse is inhibited, 
which sets an upper limit to the degree of rotational ex¬ 
citation with periodic pulse trains. Instead of increasing 
linearly with N, the molecular angular momentum os¬ 
cillates between J = 0 and the centrifugal limit. The 
oscillatory behavior has been described theoretically as 
an analogue of Bloch oscillations in solids [2^ and re¬ 
cently demonstrated experimentally in a gas of rotating 
molecules m exposed to a pulse train of a moderate 
strength, A = 8, Pn < 40. 

Here, we investigate the rotational excitation of oxygen 
by a pulse train of much higher cumulative kick strength, 
A = 20, Pat < 140. A rigid O 2 molecule would have 
been excited to the quantum states with J » 140, higher 
than that produced by an optical centrifuge (J « 120 
[21]). Yet, owing to the centrifugal limit, our high-P/v 
pulse train fails to populate rotational levels with J > 
17. We use state-resolved coherent Raman detection to 
explicitly demonstrate the effect of centrifugal distortion 
as the main reason for the inhibited rotational ladder 
climbing. 

From the analysis of the rotational Raman spectrum 
as a function of the timing of pulses in the train, one 
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FIG. 1: (color online) Experimental setup for the generation 
of long sequences of amplified femtosecond pulses. See text 
for details. 

gains an insight into the possibility of depositing larger 
amounts of angular momentum despite the centrifugal 
limit. To achieve this goal, we exploit fractional quan¬ 
tum revivals [22] and modify the pulse sequence so as to 
expose the molecules to four rotational kicks per T^ev 
Such an optimized train allows us to “pack” more pulses 
within the limited amount of time before the centrifugal 
distortion suppresses further excitation. By fine tuning 
the timing around each fractional and full revival, we ex¬ 
tend our reach from J « 17 to J « 29, utilizing more 
efficiently the cumulative strength of a long train with 
28 pulses. 

Driving molecules to the highly coherent state in a 
dense medium makes the latter a strong light modula¬ 
tor (for recent reviews, see [23|[21]). In this work, the 
effect is expressed through the generation of frequency 
sidebands via the cascaded rotational Raman scattering 
[25] . The molecular phase modulation (MPM) imparted 
on a probe pulse is maximized when the probe delay co¬ 
incides with a full revival of the rotational wave packet, 
as well as if its timing is close to that of a fractional re¬ 
vival |26j . In the case of a long pulse train optimized 
for the efficient rotational excitation, the spectral band¬ 
width is increasingly broadened from pulse to pulse, as 
has been theoretically studied in m and experimentally 
demonstrated here. 


II. SETUP 

Our experimental setup for the generation of long pulse 
sequences is shown in Fig. (technical details are dis¬ 
cussed elsewhere [2H])- Briefly, we use a Ti:Sapphire 
femtosecond laser system (SpitFire Pro, Spectra-Physics) 
producing pulses of 120 fs full width at half maximum 
(FWHM) at a central wavelength of 800 nm, 1 KHz rep¬ 
etition rate and 2 mJ per pulse. The pulses are first sent 
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FIG. 2: (color online) Temporal profile of a long periodic pulse 
train (a) and a dense aperiodic pulse train (b), 20 pulses each 
(note different time scale). The four shades of red correspond 
to one of the four different pathways through the two nested 
Michelson interferometers. 


through a pulse shaper, built in a standard ‘4/’ geom¬ 
etry [22] ■ Its spectral resolution of 0.04 nm dictates an 
upper duration limit for the shaped pulse train around 
50 ps. The shaper is followed by two Michelson interfer¬ 
ometers generating four copies of the shaped train, which 
are distinguished by four different colors in Fig. [^ 

Producing a long sequence of pulses, whether with a 
pulse shaper or an interferometer, is typically accompa¬ 
nied by a significant sacrifice of the total energy in the 
final pulse train. To compensate these energy losses, we 
amplify the train with a multi-pass amplifier (MPA), as 
shown in Fig.[^ The shaped pulses are passed four times 
through a Ti:Sapphire crystal, pumped by a Neodymium 
YAG laser (Powerlite Precision II, Continuum, 800 mJ at 
532 nm and 10 Hz). To avoid optical damage and non¬ 
linear propagation effects in air, we shape and amplify 
chirped pulses before compressing them with a grating 
compressor in front of the gas cell. The combination of a 
high resolution pulse shaper, two Michelson interferome¬ 
ters and a MPA allows us to create an arbitrary number 
of pulses within a time window of about 250 ps, with each 
pulse reaching energies of 100 ^J. As a result of the non¬ 
linear amplification process, the pulses attain a shorter 
duration around 100 fs. The standard deviation in the 
pulse amplitudes within a train is about 20%. 

We characterize the produced pulse train by cross- 
correlating it with a single transform-limited 120 fs pulse, 
which is sent through a separate pulse shaper for disper¬ 
sion compensation (Fig. [^. Two examples of such cross¬ 
correlation are shown in Fig. [^ illustrating the shaping 
capabilities used in this work. Typically, a periodic pulse 
train of five pulses with a period T4 « Trev is created by 
the pulse shaper. The interferometers are used to create 
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Rotational state J 

FIG. 3: (color online) (a) Experimental setup for coher¬ 
ent rotational Raman spectroscopy implemented in this work, 
(b) Raman spectrum of room temperature oxygen gas after 
the rotational excitation with a single weak transform-limited 
pulse. 

four copies of this train that can either be added consec¬ 
utively one after another to form a long periodic pulse 
train of 20 pulses or interleaved with variable timings Ti, 
T 2 and T 3 = Ti+ T 2 as shown in Fig. [^a) and Fig. ib), 
respectively. The latter approach enabled us to generate 
dense aperiodic pulse trains optimized for the maximum 
rotational excitation, as described below. 

A beta barium borate (BBO) crystal is used to produce 
probe pulses with the center wavelength of 400.8 nm. 
Probe pulses are linearly polarized at 45° with respect to 
the pulses in the train and are combined with the latter 
on a dichroic beamsplitter. All pulses are focused into 
a cell containing room temperature oxygen at variable 
pressure, as shown in Fig. [^a). Molecular alignment, 
produced by strong pulses in the train, results in the op¬ 
tical birefringence of the gas and a corresponding change 
in the probe polarization [30]. The latter is detected by 
passing the output probe light through an analyzer set 
at 90° with respect to the initial probe polarization. An 
optical translation stage is used to scan the probe delay. 

We use coherent Raman scattering to detect the ro¬ 
tational coherence, induced by a pulse train in oxygen 
molecules. When a weak probe pulse follows the train 
through the cloud of rotating molecules, its spectrum 
acquires Raman sidebands shifted up or down with re¬ 
spect to the central frequency. The shift of each Ra¬ 
man peak, Awj, is equal to the frequency spacing be¬ 
tween the rotational levels |J) and |J -I- 2), while its 
magnitude is proportional to the square of the rotational 
coherence pyj +2 [S]. Using a spectrally narrow probe 
(AA < 0.15nm FWHM) enables us to resolve individual 
rotational states of oxygen and to determine the shape 
of the excited rotational wave packet. Figure |^b) dis¬ 
plays such a Raman spectrum taken after the rotational 
excitation by a single weak kick. In this case, the spec¬ 
trum reflects the Boltzmann distribution among the rota¬ 
tional states of molecules. The detected wavelength shift 
can be translated to the rotational quantum number J, 
shown along the lower horizontal axis. Due to the nuclear 


spin statistics of oxygen, J can take only odd values. All 
state-resolved Raman spectra presented in this work were 
taken with the probe pulse arriving immediately after the 
last pulse in the train. 


III. RESULTS 

First, we investigate the rotational coherences created 
by a periodic pulse sequence in a thermal ensemble of 
oxygen molecules. Figure l^a) is a 2D plot showing 
the observed Raman peaks (color coded from dark to 
bright red) after a train of five relatively weak pulses 
(2 • lO^^W/cm^ or P = 0.5) as a function of the train pe¬ 
riod. The apparent pattern of peaks can be interpreted 
using the following arguments, based on the perturbative 
regime of light-molecule interaction. The first laser kick 
in the train induces multiple Raman transitions between 
the rotational states | J) of a molecule, creating a coher¬ 
ent rotational wave packet 4'(t) = |J), 

with the complex amplitudes cj and rotational energies 
Ej. Consider a wave packet of only two states, | J) and 
IJ -I- 2). It can be assigned a classical rotation period 
Tj = 2h (P ,/+2 — Ej)~^ with Planck’s constant h. Given 
the symmetry of an oxygen molecule, the wave function 
of such a coherent superposition repeats itself every inte¬ 
ger multiple of half-rotations, i.e. at times Tj = Njtj/2 
{Nj G N) indicated by light blue markers in Fig. Ha)- 
If the next laser pulse arrives at this time, it “kicks” 
the molecule in the same direction as the previous pulse, 
enhancing the rotational excitation of the corresponding 
(I J), jJ -I- 2)) wave packet. 

For a rigid rotor, Ej = hcBJ{J + 1) and the rota¬ 
tional period becomes tj = T^eviJ + 3/2)“^, with the 
J-independent revival time Trev = il2,cB. Hence, the 
number of half-rotations of an excited wave packet dur¬ 
ing its free evolution between two consecutive laser kicks 
is Nj = T/Tj-e-v X {2J 4-3). At the timing known as the 
“quantum resonance”, when T = Trev, all wave packets 
perform an integer number of half-rotations with Nj = 
2J + 3. In Figurej^a), these times Tj = {2J + 3) x tj/2, 
which are labeled with light blue squares, make up an al¬ 
most vertical “trajectory” originating at the revival time 
for low values of the rotational quantum number. Thus, 
a resonant pulse train with the time period T = Trev is 
equally efficient in exciting most molecules in the ensem¬ 
ble, regardless of their angular momentum [TS]. This is 
reflected by the highest number of Raman peaks around 
Trev (11.67 ps for ^® 02 ) marked by the vertical solid line 
in Fig. [4]; a). 

Away from the quantum resonance, each rotational 
wave packet completes an integer number of half¬ 
rotations at different time moments. For instance, in¬ 
creasing (decreasing) Nj to the next integer value 2 J-l-4 
(2 J-l-2) results in a curved trajectory marked with trian¬ 
gles (circles), along which no two rotational wave packets 
corresponding to two different values of J are in phase 
simultaneously. The agreement between the calculated 
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FIG. 4: (color online) (a) State and time resolved Raman 
spectrogram of (color coded from dark to bright red) 

after the excitation by a sequence of five pulses with a period 
T scanned around the rotational revival time T^ev (vertical 
solid line). The vertical shaded band represents the length of 
our Gaussian pulses (FWHM), whereas the shaded horizonal 
band covers the rotational quantnm nnmbers (shown on the 
right vertical axis) corresponding to the thermally popnlated 
rotational states (popnlation higher than 10 %). Light blue 
markers indicate the time moments at which a coherent rota¬ 
tional wave packet consisting of two states, | J) and | J -1-2), 
completes an integer multiple of half rotations, Nj. The cen¬ 
tral trajectory (squares) corresponds to Nj = 2J -|- 3, while 
the two neighboring sets of circles and triangles represent the 
cases of Nj = 2J + 2 and Nj = 2J + A, respectively (see text 
for details), (b) Raman spectra after a periodic sequence of 
20 pulses with T = Trev - black solid; T = 0.996 • Trev - blue 
dashed; T = 1.004 • Trev - green dotted. 


times Tj and the peaks in the measured Raman spec¬ 
tra supports the provided interpretation of the coherent 
accumulation in the rotational multi-pulse excitation. 

Even when the conditions of the quantum resonance 
are met for the lower rotational states, they no longer 
hold at higher J’s owing to the centrifugal distortion 
of the molecular bond. With increasing J, the cen¬ 
trifugal term in the rotational energy Ej = hc[BJ{J -h 
1) — DJ^{J -h 1)^] (with D being the centrifugal con¬ 
stant) becomes non-negligible, making the resonance J- 
dependent and, therefore, impossible to satisfy for all 
quantum states simultaneously. In Figure |^a), this ef¬ 
fect is seen through the apparent curving of the reso¬ 
nant trajectory (light blue squares) away from the ver¬ 
tical line at T = T^ev above J « 15. As a result, the 
efficiency of the accumulative rotational excitation by a 
resonant pulse train (T = T^ev) deteriorates with growing 
J. To demonstrate the described centrifugal limit we in¬ 
creased the number of kicks as well as their strength. 


Fig. gb) shows the detected Raman signal generated 
by the periodic resonant train of 20 strong pulses with 
pulse intensity of 3 • lO^^W/cm^ per pulse (black solid 
line). Despite the increased cumulative kick strength of 
P 20 « 140, the highest excited level remains significantly 
below that value, J « 17 ^ ^ 20 - In fact, even though 
the perturbative analysis used earlier is not applicable in 
the case of strong pulses, the measured limit agrees well 
with our conclusions drawn from Figure]^ a), as well as 
with recent numerical calculations [20j and experimental 
observations [Si- 

Utilizing the perturbative picture further, one arrives 
at a simple method of extending the reach of rotational 
excitation by a periodic train of femtosecond pulses. As 
seen in Fig. |^a), by shifting the train period above the 
quantum resonance, its overlap with the time of com¬ 
plete half-rotations Tj of the rotational wave packets 
with higher J’s can be improved. However, the degree 
of such control is rather limited, because detuning the 
pulse train far from the resonance reduces the excitation 
efficiency of the initially populated low-J states. The 
limitation can be analyzed by inspecting the two shaded 
bands in Fig. |^a): the vertical one represents the tem¬ 
poral width of our pulses (FWHM), whereas the horizon¬ 
tal one covers the thermally populated rotational states. 
The length of a continuous set of blue markers under the 
vertical band tells us about the width of the created ro¬ 
tational wave packet. On the other hand, how many of 
them are covered by the intersection of the two shaded 
areas is a qualitative indicator of the number of molecules 
in the wave packet. As the train period shifts to the right 
of the quantum resonance, the former factor increases at 
the expense of the latter. 

For the length and strength of our experimental pulse 
trains, we found the detuning of -1-0.4% to result in the 
highest enhancement of the wave packet’s width, while 
not causing a significant loss of the overall excitation effi¬ 
ciency. As illustrated in Fig. Bb). by increasing the train 
period 0.4% above the quantum resonance (green dotted 
line), we indeed shift the rotational excitation to higher 
J states. Similarly, setting the train period 0.4% below 
the resonance (blue dashed line) results in the narrower 
excited wave packet with the lower “center of mass”. The 
idea of small detunings from the resonance was used in 
m for the detection of Bloch oscillations. By tuning 
below the resonance and bringing the centrifugal limit to 
lower J values, one shortens the period of Bloch oscilla¬ 
tions, enabling their observation with smaller number of 
pulses. Bloch oscillations are not visible in this work due 
to the spatial averaging over the pump beam intensity 
profile. 

We note, that due to the power broadening of the in¬ 
dividual rotational transitions, higher pulse intensities 
should allow farther detunings from the quantum reso¬ 
nance and correspondingly broader rotational wave pack¬ 
ets. Stronger pulses effectively push the centrifugal limit 
higher. In this regime of interaction (which is beyond the 
scope of this work), the finite width of the laser pulses 
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FIG. 5: (color online) (a) Raman spectrogram for the rota¬ 
tional excitation by a sequence of two identical periodic pulse 
trains, five pulses each. The period of each train is fixed 
at Trev, while the time delay Ti between the two trains is 
scanned. The integrated signal above reveals the times of 
maximum total coherence, (b) Alignment-induced birefrin¬ 
gence signal as a function of the probe delay after a weak 
transform-limited pump pulse, (c) Raman response after the 
rotational excitation by a sequence of 20 pulses: a periodic 
pulse train with T = Trev (black solid) and an optimized pulse 
train (red dashed). See text for the optimization procedure. 
The two trains are illustrated schematically in the upper right 
corner. 

sets the ultimate limit for the reach of the rotational ex¬ 
citation. Once a molecule rotates too fast for the pulses 
to act as instantaneous kicks, it does not climb the rota¬ 
tional ladder any further. 

The effect of the centrifugal distortion grows in time 
as the wave packet spreads out. The dispersion of the re¬ 
vival times with J accumulates, making every next kick 
less and less efficient M- Hence, to fully benefit from 
the large number of pulses, one needs to apply them on 
as short of a time scale as possible. Here, we utilize 
fractional revivals [22j and create a tunable non-periodic 
pulse sequence in order to increase the efficiency of rota¬ 
tional excitation. 

To start, we demonstrate that utilizing two pulses per 
every revival time can indeed enhance the excitation, if 
the timing is chosen appropriately. For that purpose, a 
train of five pulses with a period T 4 = T^ev is produced 
with the pulse shaper. Using a single Michelson inter¬ 
ferometer, the train is then split in two parts, which are 
overlapped in space with the time delay Ti between them 


(see Fig. Hb) for the definition of time intervals Ti and 
T 4 ). Similarly to the Raman spectrogram in Fig. ia), 
in Fig. I^a) we plot the observed frequency-resolved Ra¬ 
man signal as a function of the delay Ti between the two 
interleaved periodic pulse trains. As before, strong Ra¬ 
man peaks correspond to the time moments of enhanced 
rotational excitation. At some of those time moments 
marked with solid vertical lines, slightly before or after 
the (|, ^ and |) x T^ev fractional revivals, the Raman 
response is simultaneously high for the largest number of 
rotational states, in exact analogy to the previously dis¬ 
cussed response to the periodic excitation at the quan¬ 
tum resonance. We highlight this in the upper panel by 
plotting the same Raman signal integrated over wave¬ 
length. The latter signal is proportional to IPJ.J+ 2 I 
and should not be confused with the degree of transient 
molecular alignment which is plotted in Fig. ib) for 
comparison (the alignment signal was recorded with a 
120 fs probe pulse following a single kick). One can now 
interpret the optimal Ti values in Fig.j^a) as the times of 
the maximum positive derivative of the alignment signal, 
when the majority of molecules move towards the aligned 
state. Introducing a second kick at this time results in 
the enhanced rotational excitation and strong coherent 
response. In contrast, dashed vertical lines correspond to 
the maximum negative derivative of the alignment signal, 
when the majority of molecules move towards the anti¬ 
aligned state. A second kick arriving at this time lowers 
the degree of excitation as reflected by the dips in the 
integrated coherence signal in panel (a). 

The above analysis suggests a possibility of enhancing 
the rotational excitation by adding up to six pulses per 
every revival period. Our pulse shaping scheme enables 
convenient scanning of four pulses per T^ev^ near one full 
and three fractional revivals. We set the four variable de¬ 
lay times to the following optimal values taken from the 
integrated coherence signal in Figj^a): Ti = 0.242 • Trev, 
T 2 = 0.519 • Trev, ^3 = 0.761 • Trev and T 4 = 1.004 • Trev = 
11.72 ps, i.e. slightly above the full revival as explained 
earlier in the text. We note that these optimized delays 
depend on the kick strength of the pulses in the train. 
The stronger the pulses, the higher the rotational excita¬ 
tion after each pulse, the larger the centrifugal distortion, 
the bigger the required shift from every fractional revival. 

The result of the excitation by an optimized train 
is shown in Fig. [^c). By shortening the duration of 
the pulse train, while using the same number of pulses 
(fV = 20, Fig.lKb)), we extend our reach from J=17 with 
one pulse per revival (black solid curve) to J « 29 with 
four pulse per revival (dashed red curve). Note that al¬ 
though the centrifugal limit of the periodic excitation has 
been circumvented, the efficiency of the aperiodic opti¬ 
mized pulse train is still well below its accumulated kick 
strength of Pn = 140. 

At the same time, despite the limited degree of exci¬ 
tation, the induced rotational coherence in the molecu¬ 
lar ensemble is rather strong. This can be seen through 
Raman processes of higher orders, i.e. cascaded Raman 
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FIG. 6: (color online) (a) Simplified illustration of the first three scattering orders in the process of cascaded Raman scattering. 
Grey profile represents the Boltzmann distribution of rotational population (see text for details), (b) Rotational Raman spectra 
of oxygen gas excited by the optimized train of 28 pulses (see text for the optimization procedure). The four curves correspond 
to the gas pressures of 1.7 atm (lower blue), 2.4 atm (middle orange), 5.1 atm (higher red) and 6.5 atm (top black). First, second 
and higher orders of the cascaded Raman scattering are indicated with arrows. More than 235 Raman peaks are observed at 
higher pressure values, as shown in the inset. 


scattering, especially pronounced at higher gas pressure. 
In Fig.j^b), we plot the spectra of probe pulses following 
a pulse train of 28 pulses at four different pressures rang¬ 
ing from 1.7 to 6.5 atm. The timing of the pulses in the 
train has been optimized at 6.5 atm so as to achieve the 
highest integrated coherence signal at J > 17, i.e. for 
maximizing the sum X]/>i 7 IpT 4 -i- 2 | > by means of the 
same optimization procedure as described above. 

With an optimized pulse sequence, we count up to 250 
Raman peaks. Note, that this number does not reflect 
the rotational quantum number reached by the end of 
the excitation process (indeed, at such high J values, the 
separation between the consecutive peaks should have 
decreased significantly due to the centrifugal distortion). 
Rather, large frequency shifts are caused by the cascaded 
Raman scattering giving rise to multiple coherent side¬ 
bands [5S]. The process is illustrated schematically in 
Fig. j^a) . A strong pump pulse (solid orange arrow) cre¬ 
ates multiple coherences pjj +2 via several two-photon 
Raman processes (for clarity, only one Raman transition 
with an average frequency shift A is shown in the hg- 
ure). The shape of the resulting Raman spectrum, which 
consists of a series of n red-shifted peaks, reflects the 
initial thermal distribution of molecules among the rota¬ 


tional states (grey profile). In the second-order process, 
all n emitted photons (dashed red arrow) re-scatter off 
the induced coherences, giving rise to the secondary set of 
'n? Raman lines (dotted brown arrow), centered around 
2A from the input laser frequency. The third-order pro¬ 
cess results in yet another red shift by A (dotted brown 
to dash-dotted black) and so on. This picture explains 
the repetitive broad spectral features in Fig. |^b), e.g. 
those marked as “higher orders”, resembling the Raman 
spectrum observed after a single weak pulse in Fig. ib). 
The features are separated by the frequency shift A cor¬ 
responding to the peak of the Boltzmann distribution. 
The number of Raman processes in each order m of the 
cascaded scattering grows as n™. We note, that the out¬ 
put spectrum consists of uniformly spaced Raman lines, 
owing to the repetitive scattering off the same set of co¬ 
herences, induced by the strongest hrst-order interaction. 

We verified experimentally that the conversion effi¬ 
ciency into the higher-order sidebands increases with the 
intensity of the pulse train and the number of molecules 
in the interaction region, both leading to the stronger 
rotational coherence in the system. Raising the density 
of molecules by increasing the gas pressure revealed a 
Bessel-like dependence of the strength of the Raman side- 
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FIG. 7: (color online) (a) Time-resolved Raman spectrogram 
of oxygen gas excited by the optimized train of 28 pulses (see 
text for the optimization procedure) under the pressure of 
6.5 atm. (b) Raman spectra at the time moments indicated 
with arrows in the spectrogram: before the first kick (black 
solid), at the time of the first kick (blue dashed), and at the 
time of the third kick (red dotted). 


bands on the scattering order number |25j . For instance, 
as can be seen in Fig. [^b), at P = 5.1 and 6.5 atm, 
the second-order coherences exceed those induced by the 
first-order scattering process. 

Until now, we have used narrowband probe pulses 
for the state-resolved detection of the observed spectral 
broadening driven by the Raman transitions of multiple 
orders. In the time domain, the process can be described 
as the transient molecular phase modulation (MPM) ow¬ 
ing to the periodically modulated refractive index of the 
medium. If a femtosecond pulse coincides with a full or a 
fractional wave packet revival, when the phase modula¬ 
tion is maximized, its frequency bandwidth is broadened 
[26] . The broadening effect has been theoretically pre¬ 
dicted to accumulate from pulse to pulse in a periodic 
sequence of pulses, as long as the train period is equal to 


the rotational revival time m- 

To demonstrate this accumulative broadening, we scan 
the time delay of a 120 fs probe pulse across the opti¬ 
mized pulse train described earlier, while recording the 
probe spectrum. The results are plotted in Fig. [^a) as a 
function of the probe delay with respect to the first pulse 
in the train. As expected, the spectrum of the probe 
pulse remains unchanged unless the latter coincides in 
time with a full or a fractional revival of the wave packet. 
The further down the pulse train we probe, the broader 
and more red-shifted the probe spectrum becomes, as il¬ 
lustrated in Fig. mb). 

IV. SUMMARY 


We have explored the possibilities of efficient rotational 
excitation of diatomic molecules by long trains of intense 
femtosecond pulses. In the case of a periodic pulse se¬ 
quence tuned to the optimal timing of a quantum res¬ 
onance, the excitation strength is shown to be severely 
limited by the centrifugal distortion. We have demon¬ 
strated that the centrifugal limit can be shifted up to 
slightly higher values of angular momentum by properly 
choosing the train period. 

To further extend the degree of rotational excitation, 
we have developed an optimization procedure which pro¬ 
duced non-periodic pulse sequences capable of creating 
broad rotational wave packets surpassing the centrifugal 
limit. The optimization has been demonstrated with a 
train of twenty-eight strong femtosecond pulses, with ev¬ 
ery pulse precisely timed around the corresponding frac¬ 
tional revival of oxygen. Even though the rotational lev¬ 
els, accessed with this technique, are still lower than those 
expected in the absence of the centrifugal distortion, the 
induced rotational coherence is enhanced significantly 
in comparison to the single-pulse rotational excitation. 
We have demonstrated this enhancement by observing 
the spectral broadening of probe light due to the strong 
molecular phase modulation or, equivalently, high-order 
Raman scattering. The developed method may prove 
useful for creating ensembles of highly aligned molecules 
under ambient conditions. 
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